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We study the optical properties of gapped graphene in presence of a magnetic field. We consider 
a model based on the Dirac equation, with a gap introduced via a mass term, for which analytical 
expressions for the diagonal and Hall optical conductivities can be derived. We discuss the effect of 
the mass term on electron-hole symmetry and tt-tt* symmetry and its implications for the optical 
Hall conductivity. We compare these results with those obtained using a tight-binding model, in 
which the mass is modeled via a staggered potential and a magnetic field is included via a Peierls 
substitution. Considering antidot lattices as the source of the mass term, we focus on the limit 
where the mass term dominates the cyclotron energy. We find that a large gap quenches the effect 
of the magnetic field. The role of overlap between neighboring vr orbitals is investigated, and we 
find that the overlap has pronounced consequences for the optical Hall conductivity that are missed 
in the Dirac model. 
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I. INTRODUCTION 



While already the subject of a Nobel Prize in physics, 
research in graphene^i^ a single two-dimensional sheet 
of carbon first isolated in 2004)^ seems to show no signs 
of slowing down. Initial studies have focused on the 
unique electronic properties of pristine graphene;^ such 
as, e.g., its extreme electron mobility;^ and a remarkably 
large cyclotron gap, which has led to integral quantum 
Hall measurements at room temperature^ This feature 
of graphene arises due to its linear band structure near 
the Fermi energy, which leads to an unconventional half- 
integer quantum Hall effect due to the existence of a Lan- 
dau level at zero energyji"— The field has recently ma- 
tured to a point where a large part of the focus has shifted 
to the possible applications of graphene, in fields as di- 
verse as transistors,— solar cells,— hydrogen storage,— 
and touchscreen devices.— This inevitably draws light 
to one of the serious drawbacks of graphene, namely its 
semi-metallic nature. Several ways of opening a gap in 
graphene have been put forth. If graphene is cut in nar- 
row ribbons, so-called graphene nanoribbons, quantum 
confinement effects will induce a band gap, the size of 
which depends on the width of the ribbon as well as 
the intricacies of the ribbon edgesJ^"— Adsorption on 
graphene of hydrogen has proven to be a very efficient 
way of inducing substantial band gaps in graphene, with 
fully hydrogenated graphene, termed graphane, exhibit- 
ing a band gap of several electron volts J^ii^ We have pre- 
viously suggested another means of rendering graphene 
semiconducting, by creating a periodic array of antidots, 
termed graphene antidot lattices<^i2i The source of the 
antidots may be actual perforations of the graphene sheet 
as envisioned in the original proposal, or it may be, e.g., 
patterned hydrogen adsorption.— Transistors based on 
graphene antidot lattice have successfully been fabricated 
and characterizedi^i^ 

Irrespective of the specific mechanism responsible for 
the band gap, the simplest and most general way of 



modeling it is through a mass term. Carriers near the 
Fermi level of graphene can be described via a Hamil- 
tonian that resembles the Dirac Hamiltonian of massless 
neutrinos r^ which has led to the term massless Dirac 
fermions being applied to the low-energy excitations of 
pristine graphene. Adding a term that acts oppositely on 
each sublattice is equivalent to adding a mass to these 
otherwise massless quasiparticlcs and consequently in- 
duces a band gap of twice the magnitude of the mass 
term. This results in so-called gapped graphene. 

While the focus of much research is now on inducing 
a band gap in graphene, the field remains at a point 
where the fundamental properties of the gapped struc- 
tures remain to be investigated. With this in mind, we 
focus in this paper on investigating the magneto-optical 
properties of gapped graphene. Previous studies of such 
properties have considered excitonic effects as the source 
of the mass term^i^ we have in mind structures such 
as graphene antidot lattices, for which substantial band 
gaps can be induced. Focus will thus largely be on the 
limit in which the mass term dominates or is compara- 
ble to the cyclotron energy. While the Dirac equation 
provides a reasonable approximation of the low-energy 
structure of graphene, we extend previous studies of 
magneto-optical properties by comparing with results of 
tight-binding models. Of particular interest is the effect 
of including the overlap between neighboring n orbitals, 
which is of the order s = 0.1 and thus hardly insignifi- 
cant. It is well-known that the overlap breaks electron- 
hole symmetry in the spectrum of graphene^ which is of 
significant consequence for the magneto-optical proper- 
ties of gapped graphene. In particular, we demonstrate 
that the broken symmetry results in a much richer struc- 
ture of the optical Hall conductivity of gapped graphene. 
Also, a non-zero optical Hall conductivity is found, even 
for a chemical potential sitting in the middle of the gap, 
where the electron-hole symmetry inherent in the other 
models dictates a vanishing optical Hall conductivity. 

The article is organized as follows. In section [H] we 



first present the Dirac model, and discuss some intriguing 
features of the Dirac model related to valley asymmetry 
before presenting analytical expressions for the optical 
conductivities. These expressions were first derived in a 
slightly different manner by Gusynin et al. in Ref. 1271 . 
We then present the tight-binding model, emphasizing 
the role of overlap between neighboring n orbitals. We 
generalize previous results regarding the role of overlap 
in pristine graphene to the case of gapped graphene. In 
section lllll we present the results of the different methods. 
We first focus on the diagonal optical conductivity and 
then move on the optical Hall conductivity. We compare 
the three different methods and discuss the effect of the 
mass term, particularly in the regime where the mass 
term dominates the cyclotron energy. We illustrate the 
particular importance of including overlap in the tight- 
binding model in relation to the optical Hall conductivity. 
Finally, in IIVI we summarize our findings. 



II. MODELS 

We consider gapped graphene, i.e. including a mass 
term, and include a magnetic field applied perpendicular 
to the graphene plane, B = Bz. We choose the Landau 
gauge, A = Bxy, and assume B > throughout the 
article. To calculate the optical response, we use the 
Kubo-Greenwood formula 
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where the sum is over all states, each described by a 
complete set of quantum numbers, Vn. Here, E^^^i,^ = 
Evm ~ ^I'n 1 while n{E) is the Fermi distribution func- 
tion, A is the unit cell area, hV represents a broadening 
term, and we have normalized with the zero-frequency 
graphene conductivity ctq = e^ /{AK). Also, we have in- 
troduced H"^ j^^ as the a-component of the canonical 
momentum matrix elements 11^^^^^ — {vn\a\vm), with 
n = p + eA. To evaluate the momentum matrix el- 
ements we will make use of the commutator relation 
li = imh-^[H,Y\. 

In order to make the comparison with tight-binding 
results more clear, and to illustrate in detail some of the 
intriguing properties of the model, we will proceed by 
deriving analytical results of the optical conductivity of 
gapped graphene, in a model based on the Dirac equa- 
tion. This will serve to highlight the differences and sim- 
ilarities between the Dirac and the tight-binding treat- 
ment of the problem, and provide us with some analyti- 
cal expressions with which to compare the tight-binding 
results. We stress that the analytical expressions for the 
optical conductivity have already been derived previously 
by Gusynin et al. in Ref. [23, although with a slightly 
different approach than the one we will adopt. We thus 



repeat some of these results as well as others by Jiang 
et al^, in order to assist in the comparison with the 
tight-binding results later. Also, the derivation will serve 
to clarify the effect of the mass term on the symmetry 
between electron and hole states, which has significant 
implications for the optical Hall conductivity. 



A. Dirac equation 

The low-energy excitations of graphene near the K val- 
ley are well described by a Dirac Hamiltonian^^ Hk = 
Wi?n • <T -|- Acr^, with (T = (axyCFy) where cr^ are Pauli 
spin matrices. Here, we have included a mass term Acr^, 
which breaks the sublatticc symmetry and adds a band 
gap to graphene of size 2A. The two incquivalent val- 
leys of graphene are related via time-reversal symmetry 
and a Hamiltonian valid near the K' valley is thus ob- 
tained by substituting cr — > cr*. In matrix notation, the 
Hamiltonian thus reads 
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Here, the upper (lower) sign corresponds to the K {K') 
valley, and we have introduced H± = vp{pxHpyHeBx), 
with vp the Fermi velocity of graphene. The eigenvalues 
of this Hamiltonian are^^ 



E'^ = sgn{n)^A^ + h^uM~SnoA, n = 0,±l, 
sgn(n)^A2 + h^oj^\n\ + (5„oA, n = 0,±l. 
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where ujc = v2wf/^s, with the magnetic length Ib ~ 
■\/ft/(ei?), and we define sgn(O) — 0. This results suggest 
that in the presence of a mass term, particle- hole symme- 
try is no longer retained in each individual valley, where it 
is broken by the single Eq = ±A eigenstatci^ In addition 
to the absence of particle-hole symmetry in the energy 
spectrum, we find that the general tt — tt* symmetry be- 
tween electron- hole pairs of eigenstates is broken, because 
the mass term breaks sublattice symmetry. Letting AJ^^ 
and B^j, denote the individual spinor components of the 
eigenstates, we find that while in the absence of a mass 
term, the tt—tt* symmetry relates a given hole state to its 
corresponding electron state via ^^„ i^ = {—^nk^E^i.)^, 
the mass term results in a relation instead between states 
in opposite valleys, ^5„ ^ = {B^k^-^nk)^- We will see 
that this has particular consequences for the optical Hall 
conductivity in gapped graphene. 

To evaluate the momentum matrix elements we note 
that the commutator relation yields 11 = mvpty while 
n^' = (n^)*. This yields the transition rule \m\ = 
|n| ± 1 for optical transitions \n) — !• \m). We find that 
the optical conductivity tensor can then be written in the 
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FIG. 1: (Color online) Transitions contributing to the optical 
Hall conductivity in gapped graphene, in the case of zero tem- 
perature and chemical potential. The weights of each transi- 
tion is indicated, with color indicating the sign. Note that for 
vanishing gap, Ak = B^ and the Hall conductivity vanishes 
within each valley. This does not hold in the gapped case, 
where instead contributions from opposite valleys cancel. 



form 
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where all energies arc now in units of the hopping element 
t, 17 = fta;/|t| and 7 = /ir/|t|, and we have introduced the 
relative magnetic flux (j> = Ba? /^q, with the magnetic 
flux quantum <I>o = Trh/e. Also, we define 
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where ef^ — E'^/\t\, with k G {K,K'}. Also, we have 
introduced a,'; = sgn{E^ + A)^{1 + A/E;^)/2 and b^ = 
i\/{l — A/E.^)/2. In terms of the relative magnetic flux 
e« = sgn(n)^(52 + (37r/2)0|n| =F '^no'^- Note that the ex- 
pressions in Eq. ([5]) only differ for the different valleys if 
the n — state is involved in the transition. 



The expression Eq. ((4]) for the optical conductivity 
highlights an interesting feature of gapped graphene. Or- 
dinarily, at zero temperature and chemical potential, 
TT — TT* symmetry means that contributions from tran- 
sitions |— n)— >■ \m) and |— m)— > \n) exactly cancel in the 
sum for the off-diagonal optical conductivity.— Letting 
A = we indeed find |a„p = |6„p and thus an exact 
cancellation of terms in the sum in Eq. ((3]), even if re- 
stricted to just a single valley. However, for A 7^ this 
equality no longer holds and the lack of tt — tt* symme- 
try results in a non-zero optical Hall conductivity within 
each individual valley, as discussed also in the DC case 
by Jiang et al^ Instead, the modified symmetry induced 
by the mass term between electron and hole states in 
opposite valleys means that contributions from transi- 
tions |— n) — > I to) in one valley are canceled by transi- 
tions \—m) — > \n) in the opposite valley, resulting in the 
expected axyii^) = at zero temperature and chemical 
potential, as illustrated in Fig. [T] It it thus crucial to 
take into account both valleys, and to treat the asymme- 
try of the valleys with respect to energy spectrum and 
eigenstates properly. 



In the general case, summing contributions from both 
valleys, the diagonal optical conductivity can be written 
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while the optical Hall conductivity reads 
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where 6 = A/t and we have defined eg = 6. These results are in agreement with those of Gusynin et al. in Ref. 1271 . 



while derived in a different manner. We state them again 
here for completeness. 



B. Tight-binding 

Graphene is relatively well described in a nearest- 
neighbor tight-binding (TB) framework, considering only 
TT-electrons. While commonly 7r-orbitals lp{y — R) cen- 
tered on different sites R are assumed orthogonal, we 
will consider also the case where the overlap between 
neighboring orbitals is included via the overlap integral 
s = 0.1. This turns out to have a significant impact 
on the off-diagonal term of the optical conductivity ten- 
sor. To include the effects of a magnetic field, we use 
the Peierls substitution, i.e. we adopt a minimum cou- 
pling substitution p — )■ p + eA and choose the basis set 
(r|R) = e*'^(^'""^(/3(r — R), where the phase factor is given 
as a line integral^! 
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In this basis, the Hamiltonian matrix elements can then 
be approximated as 



(Rli/IR') = (R|i/o|R'>e'^'^'^'' 



(9) 



where Hq denotes the Hamiltonian in absence of a mag- 
netic field. In a coordinate system where the x axis is 
aligned with carbon bonds, see Fig. [2l the phase factor 
becomes 

</.(R,R') = 5|-(.T + x')(zy-y). (10) 

Because of the {x 4- x') term, we are forced to choose a 
unit cell substantially larger than the Wigner-Seitz cell of 
graphene in order to retain periodicity in the problem. In 
particular, we use a rectangular unit cell of area L^ x Ly, 
with ij. = 4<I>o/(^a) and Ly = a, where a = 2.46 A is 
the graphene lattice constant. The unit cell is illustrated 
in Fig. [2l Letting N denote the number of atoms in 
the magnetic unit cell, we have A^ = 16^Q/{^/3a^) x 
B-^ ~ 316 10^ T X B-^, illustrating the disadvantage of 
the tight-binding approach, namely that very large unit 
cells are required in order to simulate realistically small 
magnetic field strengths. 

Analogously to the case of the Dirac equation, a gap 
can be introduced in the tight-binding description by 
adding a staggered potential with opposite sign on each 
sublattice, i.e. taking diagonal elements (R|i/|R) = 
Ett + (— 1)"A, with n = 0,1 for orbitals sitting on the 
A and B sublattice, respectively. Here, for generality, we 
have included the on-site energy e^r of the 7r-orbitals. 

If we take into account the overlap between neighbor- 
ing TT-orbitals, the problem becomes that of a generalized 
eigenvalue problem. However, advantage can be taken of 
the similar form of the Hamiltonian H and the overlap 
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FIG. 2: (Color online) Unit cell used in the tight-binding 
calculations. The colored circles indicate the carbon atoms 
included in the unit cell, with red and blue indicating different 
sublattices. The dashed rectangle indicates the fundamental 
unit cell, which is repeated in order to ensure translational 
symmetry. The lattice vectors of the enlarged unit cell are 
denoted A^ = L^it and Ay = Lyy. Also shown are the three 
nearest neighbor vectors Sn. 



matrix S in the chosen basis if we consider just nearest- 
neighbor coupling terms. In particular, the generalized 
eigenvalue problem can be written in the form 



(e, + A)I tF 
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where I is the identity matrix, the eigenvectors have been 
separated into components on different sublattices, and 
we have introduced the nearest-neighbor transfer integral 
t, evaluated in the absence of a magnetic field. The simi- 
lar form of the matrices allows us to rewrite the problem 
as an ordinary eigenvalue problem for either sublattice, 
and in this way arrive at an equation 



(E - e^f - A2 {Eo - e^f 



[Es - tf 
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relating the eigenvalues E to the eigenvalues Eq, obtained 
by ignoring the overlap. This is readily solved to yield 
the explicit relation 

^-^37^W^A^xi*^o-.(Co^-A^) 

+ {e.-el){t'-s^^l-A^)t''], (13) 

where ^q = Eq — e^. Here, the sign of the square root 
should follow the sign of the energy Eq. We can thus 
immediately determine the eigenvalues with the over- 
lap included, without solving an additional generalized 
eigenvalue problem. Furthermore, we use the same eigen- 
vectors, only properly orthonormalized according to the 
overlap. While using Eq. ((T3)) to calculate the eigenval- 
ues with overlap is exact, this way of determining the 
eigenvectors constitutes an approximation. However, we 
have confirmed numerically that this has no discernible 
effect of the calculated optical conductivities so long as 
sA< t. 

In addition to offering a decrease in computation time, 
the similarities of the H and S matrices allows us to 



greatly simplify the calculation of the momentum m.atrix 
elements. Following ideas similar to those used by Sandu 
in Ref. |32| for evaluating the momentum matrix elements 
of ordinary graphene, we consider the Hamiltonian of the 
ordinary eigenvalue problem H5 = S^^H. The inverse 
of the overlap matrix can be written 
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(e^ + A) I - stFFt (i-s(e^-A))F 
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The form of this matrix makes it evident how the inclu- 
sion of overlap effectively results in next-nearest neighbor 
coupling. In the crystal momentum representation, the 
momentum operator is proportional to Vki?, which for 
the Hamiltonian above becomes 



VkHs = 
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If we now assume sA <^ t and let e-,^ = Q this matrix 
simplifies to 



VkHs = t 
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(17) 
For gapped graphene in the absence of a magnetic field, 
the second term is proportional to the unit matrix and 
thus does not contribute to the momentum matrix el- 
ements, meaning that the momentum operator is unal- 
tered by the inclusion of overlap in the model, as is the 
case for ordinary graphene with no mass terni;^ How- 
ever, in the present case, where a magnetic field is in- 
cluded, such a simplification can no longer be made. In- 
stead, we note that the second term relies on two small 
effects, that of the overlap and that of the magnetic field. 
Consequently, we ignore the effect of the magnetic field 
on the matrix elements of F in this term, which allows 
us to disregard the second term entirely, as it no longer 
contributes to the momentum matrix elements. Effec- 
tively, we are thus ignoring non-orthogonality of our ba- 
sis set and calculating the momentum matrix elements 

via n„„ = i{m/h)Hmn{'R-ra - Rn)- 

In the numerical implementation of the tight-binding 
model we employ a simple equidistant fc-point mesh. Our 
choice of unit cell results in negligible dispersion along 
kx , so we only discretize along the ky direction, ensuring 
that the folded high-symmetry points are included in the 
discretization. We have verified that including more than 
a single k^ point has no influence on the results for any 
realistic values of the magnetic field. 



III. RESULTS AND DISCUSSION 

We now turn to the results for the optical conductiv- 
ity. Unless otherwise stated, from hereon we will use 
parameters Ett = 0, s = 0.1 and t = —2.2 eV, we in- 
clude a broadening term KT = 0.05 eV and fix the tem- 
perature at 300 K. Our choice of the value of the hop- 
ping term t is motivated by the transition energy of the 
saddle point resonance in ordinary graphene, which oc- 
curs at 2\t\ = 4.4 eV^i^ We note that this value of the 
hopping term results in a Fermi velocity that does not 
match experimental results. Without introducing inter- 
action terms in the form of excitonic effects, one cannot 
account both for the location of the saddle point reso- 
nance and the proper value of the Fermi velocity. In the 
end, the exact choice of hopping term does not have a 
qualitative infiuence on the results we arrive at. 

Because of the inherent problems of the Peierls substi- 
tution for the tight binding method, namely the inconve- 
nient scaling of the unit cell with magnetic field, we are 
limited to quite large magnetic fields. Therefore, most 
results will be for a substantial magnetic field of 78 T. 
However, we expect the qualitative features and the con- 
clusions to be valid at more realistic magnetic fields. This 
claim will be substantiated when we compare the tight 
binding results to those obtained using the Dirac equa- 
tion. 



A. Diagonal conductivity 

We now consider gapped graphene with a mass term 
A = 0.2 eV. We fix the chemical potential at midgap, 
/i = 0, and use T = 300 K. In Fig. [3] we show the resulting 
diagonal optical conductivity calculated using DE and 
TB with and without overlap. The results from all three 
methods indicate a very clear absorption edge, in excel- 
lent agreement with the DE solution, which predicts the 

first resonance peak at huiQ = A ( 1 -I- a/1 -|- /i^w^/A^ j . 

Focusing attention on photon energies below 1.5 eV, we 
find almost perfect agreement between all three methods, 
as illustrated in the inset of Fig. [31 Here, the resonances 
fall perfectly on the transition energies predicted by the 
DE model, as indicated by the vertical lines in the inset. 

However, a notable feature of the full range of energies 
is the peak near 2\t\ = 4.4 eV, which stems from the van 
Hove singularity at the M point of graphene. Naturally, 
such a band structure feature is not reproduced in a DE 
model. This M-point resonance is unaffected by the ad- 
dition of a mass term, so long as A ^ t. We note that 
the overlap results in a slight blueshift of the resonance 
peak, so that it occurs at 2|i|/(l — s^). 

In Fig. 21 we show the diagonal optical conductivity 
in the case where the chemical potential sits on top of 
the lowest Landau level at /x = A. The DE model sug- 
gest that in this case the lowest resonance, fejo, will split 

into two resonances at huj± = A ( ^Jl + h?uj^/A'^ ± 1 j— . 
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FIG. 3: (Color online) Diagonal conductivity a^x in units of 
the DC graphene conductivity ctq, as a function of photon en- 
ergy fkd. Results are shown for a magnetic field B = 78 T. 
Thick lines indicate results from the TB model including over- 
lap, while thin lines are without overlap. Dashed lines are DE 
results. Blue (red) coloring indicates the real (imaginary) part 
of the conductivity. Note the resonance around hu — 4.4 eV, 
which arrises due to the van Hove signularity at the M point 
of graphene, and is thus absent in the DE results. The inset 
offers a closer view of the low-energy oscillations, illustrating 
the excellent agreement between all three methods in this en- 
ergy range. The vertical lines mark the transition energies 
predicted by the DE model. 



FIG. 4: (Color online) Diagonal conductivity a^x in units of 
the DC graphene conductivity ao, as a function of photon en- 
ergy hu. Results are shown for a magnetic field _B = 78 T. 
The chemical potential is set at the value of the mass term, 
fi — A. Thick lines indicate results from the TB model in- 
cluding overlap, while thin lines are without overlap. Dashed 
lines are DE results. Blue (red) coloring indicates the real 
(imaginary) part of the conductivity. The inset offers a closer 
view of the low-energy oscillations, illustrating the excellent 
agreement between all three methods in this energy range. 
The vertical lines mark the transition energies predicted by 
the DE model. 



0, huJcC^O.23 eV 



Other than this splitting of the lowest resonance, the 
results arc nearly identical to those obtained with the 
chemical potential at midgap. The inset in Fig. U] again 
shows excellent agreement between all three methods. 

In Fig. [S] we illustrate the effect of the mass term on 
the diagonal conductivity. We show results from the TB 
model with overlap, and with a value of the mass term 
increasing from 0.1 eV to 1 eV. We fix the chemical poten- 
tial at midgap. Considering the eigencnergics obtained 
using the Dirac equation. En = A(-\/l + {h?u}'^/A'^)\n\ ± 
1), we expect that the effect of the mass term will be to 
compress the Landau level spacing, resulting in a nearly 
continuous spectrum for A 3> hiOc- This trend is evident 
in the results of Fig. [SJ where we note that the oscilla- 
tions due to Landau level spacing are almost completely 
washed out for A = 1 eV. Also shown in the figure is the 
optical conductivity in the absence of a magnetic field.— 
Comparing these to the results with magnetic field, it is 
clear that for a sufficiently large mass term, A ^ huc, 
the results are nearly identical, indicating that a large 
mass term completely cancels the effect of the magnetic 
field on the diagonal optical conductivity. We stress that 
while the cyclotron energy in this case is so large that 
a mass term of the order of 1 eV is needed to counter 
the effect of the magnetic field, more realistic values of 
the magnetic field strength would of course lower this 
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FIG. 5: (Color online) Diagonal conductivity a^x in units of 
the DC graphene conductivity ao, as a function of photon en- 
ergy huj. Results are shown for a magnetic field _B = 78 T 
and for four different values of the mass term, all calculated 
using the TB model with overlap. Full (dashed) lines indicate 
the real (imaginary) part of the conductivity. For compari- 
son, the thin lines show the optical conductivity in absence 
of a magnetic field. The vertical lines indicate the resonance 
energies predicted by the DE model. 



threshold significantly. 

In general, we find that even at the very large magnetic 
field of 78 T used for these calculations, and at chem- 
ical potentials well above the value of the mass term, 
the results for the diagonal conductivity are in excellent 
agreement between the three methods, so long as photon 
energies are well below the Af-point resonance. We see 
no reason why such an agreement should break down at 
lower field strengths, provided that broadening docs not 
smear out the resonances, and we thus argue that the 
results we obtain will remain qualitatively the same also 
at more realistic magnetic fields, with resonances simply 
red-shifted in accordance with the lower cyclotron energy. 



B. Hall conductivity 

We now turn to the off-diagonal Hall conductivity. 
Before discussing the optical Hall conductivity, we fo- 
cus on the DC Hall conductivity a^ = <Jxy{0) as func- 
tion of the chemical potential. An intriguing conse- 
quence of the linear dispersion relation of graphene is 
the emergence of an unconventional quantum Hall effect, 
with the DC Hall conductivity quantized according to 
a'^ = — (4e^//i)(n -f- 1/2), with n a positive integer or 
zero. This unique feature has its origin in the n = 
Landau level, the degeneracy of which is half that of the 
higher Landau levels.—"^ This carries over to the case of 
gapped graphene, if contributions from both valleys are 
taken into accounti^ 

In Fig. ini we show the DC Hall conductivity as a func- 
tion of chemical potential, calculated using the three 
models. In these calculations, the temperature has been 
set to IK and the broadening term HT has been omitted, 
in order to properly resolve the quantum Hall plateaus. 
In panel (a) of the figure, we show the full range of chem- 
ical potentials considered. While there is quite good 
agreement between all three models for low values of the 
chemical potential, ^ < 0.5 eV, the breakdown of the 
Dirac model is apparent at higher chemical potentials, for 
which the non-linear part of the band structure is probed. 
This results in discrepancies in the energies of the Landau 
levels, which is apparent in Fig.[6]3. While the TB model 
without overlap obviously results in a slightly different 
Landau level structure, the values of the DC Hall con- 
ductivity still fall at the same plateaus as predicted by 
the DE model. The inclusion of overlap alters the results 
quite significantly, with any but the lowest Hall plateaus 
falling at non-equidistant values. We stress that these re- 
sults arc of course calculated for a substantial magnetic 
field, and we expect much better agreement between the 
three models when the Landau level spacing is smaller. 
The TB results show an abrupt change in behavior near 
a chemical potential of /i = (1 + s) |i| , where a°y suddenly 
changes sign. This value of the chemical potential coin- 
cides with the maximum of the band at the M point of 
gapped graphene without a magnetic field, as illustrated 
in the inset of Fig. [51 The degeneracy of the Landau 
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FIG. 6: (Color online) DC Hall conductivity in gapped 
graphene as a function of the chemical potential. The thick, 
blue lines are results of the TB model including overlap, while 
thin, red lines are without overlap. Dashed, black lines are 
the DE results. Note that in contrast to the other figures, here 
the conductivity is given in units of 4e^//i. (a) DC Hall con- 
ductivity in the full range of chemical potentials considered, 
illustrating the discrepancies between the three models and 
the abrupt change in behavior near n — {1 + s)\t\ for the TB 
model results. The inset shows the band structure of gapped 
graphene without magnetic field, (b) and (c) offer closer views 
of the regimes indicated by black rectangles in panel a. Hor- 
izontal lines in (b) indicate a'^y = — (4e^//i)(n -|- 1/2), while 
those in (c) indicate a^^y — (2e /h)n. 



levels is abruptly altered at this point, with the four- 
fold degeneracy for Landau levels below E ~ (1 + s)\t\ 
changed to a two-fold degeneracy above this energy, as 
we have confirmed by inspection of the eigenvalues of 
both TB models. This results in a new structure of the 
quantized Hall plateaus, as shown in Fig. [6t, with the 
reduced degeneracy manifesting as Hall plateaus at val- 
ues (J^y = {2e'^/h)n. As they have their origin in the 
non-linear part of the band structure, these features are 
of course entirely absent in the DE results. 

We now turn to the optical Hall conductivity. As de- 
scribed previously and illustrated in Fig.[TJ the DE model 
predicts an optical Hall conductivity identically zero at 
zero temperature and for a chemical potential fixed at 
midgap, due to the complete cancellation of conjugated 
transitions in opposite valleys. Below, we will show that 
this does not hold for the TB model with overlap, but for 
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FIG. 7: (Color online) Optical Hall conductivity a^y in units 
of the DC graphene conductivity ao, as a function of photon 
energy hoj. Results are shown for a magnetic field _B = 78 T. 
The chemical potential is set at the value of the mass term, 
^ = A. Thick lines indicate results from the TB model in- 
cluding overlap, while thin lines are without overlap. Dashed 
lines are DE results. Blue (red) coloring indicates the real 
(imaginary) part of the conductivity. Vertical lines indicate 
the resonance energies predicted by the DE model. The inset 
illustrates the transitions contributing to the Hall conductiv- 
ity. See text for more details. 
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FIG. 8: (Color online) Optical Hall conductivity a^y in units 
of the DC graphene conductivity ctq , as a function of photon 
energy /kj. Results are shown for a magnetic field i? = 78 T 
and for four different values of the mass term, all calculated 
using the TB model with overlap. The chemical potential is 
in each case fixed at the lowest Landau level, p = A. Full 
(dashed) lines indicate the real (imaginary) part of the con- 
ductivity. The inset shows a closer view of the first resonance. 
The results of the TB model without overlap is shown by 
thin, black lines, illustrating the slight discrepancies between 
the models in this regime. The DE results are practically 
identical to those of the TB model without overlap. 



now we focus on the case where the chemical potential is 
located at the lowest Landau level, /x = A. In this case, 
the symmetry argument obviously breaks down due to 
the different Fermi distribution functions of the conju- 
gated transitions, and all three models predict non-zero 
Hall conductivities. 

In Fig. [7] we show the optical Hall conductivity of 
gapped graphene with a mass term of A = 0.2 eV and the 
chemical potential at /i = A. Common to the results of 
all three models are the two pronounced resonances, that 

the DE model predict at hu;± = A (y/TThFuJlJ^ ± l) , 

similar to the lowest resonances of the corresponding di- 
agonal conductivity shown in Fig. |4l However, contrary 
to the diagonal conductivity, higher-energy resonances 
are strongly suppressed in the case of the TB model with 
overlap, and are entirely absent in the two other models. 
This can be explained quite straightforwardly by con- 
sidering the conjugate transitions and the effect of the 
overlap on the band structure. We illustrate this in the 
inset of Fig. [T] Here, the thin, red arrows indicate the 
transitions contributing to the optical Hall conductivity 
in the case of the DE model and the TB model without 
overlap. Because of perfect electron-hole symmetry in 
both these models (when including both valleys), contri- 
butions from all other transitions are canceled by their 
conjugates, the only transitions contributing being those 
involving the state at £' = A. When overlap is included 
(indicated by thick, blue arrows in the inset), electron- 



hole symmetry is broken and the conjugate transitions 
no longer cancel due to a slight difference in their transi- 
tion energies (dashed arrows in the inset). This results in 
a much richer structure of the optical Hall conductivity 
when including overlap, with additional resonances oc- 
curring due to transitions between higher-lying Landau 
levels. However, because the energy difference between 
electron and hole states induced by the overlap is only 
of the order s^ , the strengths of these resonances arc re- 
duced significantly compared to those stemming from the 
E ^ A Landau level. 

In Fig. [5] we illustrate the effect of the mass term on 
the optical Hall conductivity. We show results of the TB 
model with overlap, and increase the mass term from 
A = 0.1cVtoA = l eV, while in each case fixing 
the chemical potential at the energy of the lowest Lan- 
dau level, /i = A. As for the diagonal conductivity, we 
clearly see that the effect of the magnetic field is strongly 
reduced when the mass term dominates the cyclotron en- 
ergy. Thus, whereas at large values of the mass term the 
diagonal conductivity retains the resonance correspond- 
ing to a transition energy of roughly the size of the band 
gap, this resonance vanishes in the optical Hall conduc- 
tivity, once the mass term grows sufficiently large. In the 
inset of Fig. |8l we show a closer view of the first reso- 
nance at i?i — A, and include results of the TB model 
without overlap. While there are some small discrepan- 
cies between the two models in this regime, the overall 
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FIG. 9: (Color online) Optical Hall conductivity with a chem- 
ical potential fixed at midgap, calculated using the TB model 
with overlap. The vertical lines indicate the transition ener- 
gies predicted by the DE model. 



agreement between the models is quite good. The corre- 
sponding results of the DE model are nearly identical to 
those of the TB model without overlap. We note that the 
DC Hall conductivity is also decreased as the mass term 
increases, due to squeezing of the Landau level spacing. 
As noted above, an interesting consequence of includ- 
ing overlap in the TB model is that the electron-hole 
symmetry is broken, and exact cancellation of conjugate 
transitions at zero chemical potential no longer applies. 
We therefore expect that the inclusion of overlap in the 
TB model will result in a non-zero optical Hall conduc- 
tivity even at vanishing temperatures and the chemical 
potential fixed midgap. In Fig. [9] we illustrate this, show- 
ing the optical Hall conductivity calculated using the TB 
model with overlap, with a mass term A = 0.2 eV and 
the chemical potential at midgap. As expected, while the 
other two methods predict a zero Hall conductivity, the 
inclusion of overlap breaks the electron-hole symmetry 
and results in a non-zero optical Hall conductivity. Over- 
lap is thus a crucial ingredient for properly evaluating the 
optical Hall conductivity of graphenc. However, we note 
that because it stems from the overlap, the magnitude 
of the Hall conductivity remains significantly lower than 
the case where the chemical potential sits at the lowest 
Landau level. 



IV. SUMMARY 

We have calculated the optical conductivity tensor of 
japped graphenc in presence of a magnetic field, by em- 



ploying a Peierls substitution in a nearest-neighbor tight- 
binding model, both with and without overlap. By gener- 
alizing results from ordinary graphenc with no magnetic 
field, we have found a simple relation between energy 
eigenvalues calculated with and without overlap, allowing 
us to avoid solving a large generalized eigenvalue prob- 
lem to account for the overlap. This simplification is 
valid so long as sA <C t. We have compared the optical 
conductivities calculated using the tight-binding model 
with analytical expressions obtained in a Dirac equation 
approach)^ To elucidate the differences between these 
two approaches we have highlighted the role played by 
the mass term in redefining the symmetry between elec- 
tron and hole states in the Dirac model. This results 
in a non-zero optical Hall conductivity in individual val- 
leys even with a chemical potential fixed in the middle 
of the mass gap. However, summing contributions from 
both valleys, the symmetry induced by the mass term be- 
tween electron and hole states in opposite valleys results 
in an optical Hall conductivity identically zero. This re- 
sult carries over to the tight-binding model without over- 
lap, and in general we find excellent agreement between 
the optical Hall conductivities calculated with these two 
methods. 

Including overlap in the tight-binding model strongly 
modifies this picture, because the energy spectrum no 
longer has perfect electron-hole symmetry. We find that 
this results in a much richer structure of the optical 
Hall conductivity, while also predicting non-zero optical 
Hall conductivity even with a chemical potential fixed 
midgap. We conclude that overlap is a crucial ingredient 
for a proper, full description of the off-diagonal magneto- 
optical properties of gapped graphenc. While the optical 
Hall conductivity does show some significant discrepan- 
cies between the three models, in general we find that 
the low-energy transitions are, as expected, very well cap- 
tured in the simple Dirac model. In particular, the diago- 
nal optical conductivity at transition energies well below 
the saddle point resonance at 4.4 eV is almost identical 
between all three models. 

Finally, we find that a sufficiently large mass term, 
much larger than the cyclotron energy, effectively washes 
out the infiuence of the magnetic field on the optical prop- 
erties of gapped graphene. 
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